Synthetic LISA: Simulating Time Delay Interferometry in a Model LISA 
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We report on three numerical experiments on the implementation of Time-Delay Interferometry 
(TDI) for LISA, performed with Synthetic LISA, a C++/Python package that we developed to 
simulate the LISA science process at the level of scientific and technical requirements. Specifically, we 
study the laser-noise residuals left by first-generation TDI when the LISA armlengths have a realistic 
time dependence; we characterize the armlength-measurements accuracies that are needed to have 
effective laser-noise cancellation in both first- and second-generation TDI; and we estimate the 
quantization and telemetry bitdepth needed for the phase measurements. Synthetic LISA generates 
synthetic time series of the LISA fundamental noises, as filtered through all the TDI observables; it 
also provides a streamlined module to compute the TDI responses to gravitational waves according 
to a full model of TDI, including the motion of the LISA array and the temporal and directional 
dependence of the armlengths. We discuss the theoretical model that underlies the simulation, its 
implementation, and its use in future investigations on system characterization and data-analysis 
prototyping for LISA. 

PACS numbers: 04.80.Nn, 07.60.Ly, 95.55.Ym 



I. INTRODUCTION 



The Laser Interferometer Space Antenna (LISA) is a 
joint NASA-ESA deep-space mission aimed at detecting 
and studying gravitational radiation in the 10~ 5 -10 _1 
Hz frequency band It is expected to be launched in 
the year 2012, and to start collecting scientific data ap- 
proximately a year later, after reaching its final orbital 
configuration Q • LISA consists of three widely separated 
spacecraft flying in a triangular, almost equilateral con- 
figuration, and exchanging coherent laser beams; grav- 
itational waves (GWs) will be measured by picometer 
interferometry as modulations in the distance between 
the spacecraft. 

LISA, which will operate in a lower frequency 
band than ground-based GW interferometers, holds the 
promise of providing access to entirely new classes of 
GW sources, but it also introduces complications un- 
known to ground-based detectors, such as the complex 
signal and noise transfer functions, the problem of can- 
celing the otherwise overwhelming laser phase noise in an 
unequal-arm interferometer, the necessity of dealing si- 
multaneously with many continuous signals (including a 
confusion- noise background of galactic white-dwarf bina- 
ries), and the possibility of using multiple intcrferometric 
observables as a virtual network of GW interferometers. 
These complications hinder the analytical characteriza- 
tion of LISA's detection capabilities as a function of its 
configuration, as well as the development of data-analysis 
techniques aimed at specific GW sources. Computer sim- 
ulations will therefore play a crucial role in exploring 
LISA's performance, in obtaining insight about its op- 
timal operation, and in prototyping and testing data- 
analysis protocols. 



In this paper we report on three numerical investiga- 
tions performed with Synthetic LISA, a software pack- 
age that we developed at the Jet Propulsion Labora- 
tory to simulate the LISA science process at the interface 
of scientific and technical mission requirements. All in- 
vestigations focus on outstanding implementation issues 
for Time-Delay Interferometry (TDI), the LISA-specific 
technique currently envisaged to suppress the otherwise 
overwhelming laser phase noise by combining (with de- 
lays) the basic LISA phase measurements aboard the 
three spacecraft into composite laser-noise-free observ- 
ables (see Sec. Ill Dl for a discussion and full references). 
More in detail: 



1. In Scc. lIV Al we give the first quantitative estimate, 
based on a straight simulation, of the improvement 
in laser-noise stabilization that would eliminate the 
need for second-generation TDI for a realistic flex- 
ing LISA array 0. [!■ Q using standard Michelson 
observables. We find that an rms improvement fac- 
tor between 3 and 10 is sufficient. We give also nu- 
merical evidence of effective laser-noise subtraction 
with second-generation observables. 
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2. In Sec. lIVBl we evaluate the armlength-ranging ac- 
curacy Q that would be required for effective laser- 
noise cancellation in first- and second-generation 
TDI Michelson observables. We find that ranging 
accuracies between 30 m and 100 m (rms) are ad- 
equate when simple linear extrapolation is used to 
compute the armlengths between measurements. 

3. In Sec. IIV CI we estimate the granularity that can 
be allowed in the quantization of phase measure- 
ments while preserving effective laser-noise cancel- 
lation. Assuming white laser frequency noise ban- 
dlimitcd at 1 Hz, we find that a total of 32-34 (or 
36-38) bits per sample are needed for the Michelson 
observables of first-generation (second-generation) 
TDI. 
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We present these results as representative of the numer- 
ical experiments that become possible with state-of-the- 
art simulators such as Synthetic LISA, and we suggest 
possible directions of investigation in the final section of 
this paper. 

Synthetic LISA represents the evolution of previous 
simulation tools developed in the LISA Project 0. 
Among other improvements, Synthetic LISA is based on 
a complete model of TDI: the LISA armlengths change 
realistically with the motion of the array; the laser beams 
propagate causally; and a full set of TDI combinations 
can be generated. Synthetic LISA joins other existing 
software that simulates the LISA response to noise and 
GWs, such as the well-established LISA Simulator by 
N. Cornish and L. Rubbo Why write a new simu- 
lator, then? Being able to rely on a plurality of simula- 
tion tools allows for mutual validation and verification, 
which is crucial if implementation choices must be pred- 
icated on the results of numerical experiments. In addi- 
tion, the two simulators have a slightly different focus. 
The LISA Simulator was conceived to interface source 
simulations to data analysis, while Synthetic LISA was 
targeted to explore the interaction between LISA science 
and technology, and it must therefore operate at a lower 
level of abstraction: in particular, Synthetic LISA per- 
forms an explicit time-domain simulation of interferom- 
etry, including the cancellation of laser phase noise. On 
the other hand, it operates at a higher level of abstrac- 
tion than integrated- modeling simulations J9| : it does not 
need to model spacecraft subsystems, but rather it as- 
sumes nominal specifications of their performance. 

This paper is laid out as follows. In Sec. |n] we de- 
scribe the theoretical model of the LISA science process 
used in our simulations; in Sec. IHII we briefly discuss the 
implementation and usage of Synthetic LISA; in Sec. II VI 
we report on our main numerical experiments; and in 
Sec. E] we give our conclusions. Appendices El an d [5] 
describe, respectively, the geometric conventions and the 
treatment of noise used in Synthetic LISA. In the follow- 
ing, we set G = c = 1 unless otherwise indicated. 



II. MODELING OF A SYNTHETIC LISA 

Figure ^ is a block diagram of the LISA science pro- 
cess, as modeled in Synthetic LISA. At the top of the 
hierarchy sit the TDI observables, which represent the 
main scientific product of the mission, and which will 
be run through data-analysis algorithms to search for 
GW signals. The TDI observables are time-delayed com- 
binations of the basic interferometric measurements (y 
and z) that compare the frequencies (or phases) of the 
two lasers on each spacecraft between themselves, and 
with the lasers incoming from the other two spacecraft. 
The Doppler measurements bear the imprint of the in- 
strumental noises and of the GW signals, but the latter 
can be read off efficiently only from the TDI observables, 
which are free of the otherwise overwhelming laser phase 
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FIG. 1: A block diagram of the LISA science process. 



noise and optical-bench noise. The time-dependent ge- 
ometry of laser propagation across the LISA array influ- 
ences the effect of the LISA noises and (especially) of GW 
signals on the Doppler measurements; a precise knowl- 
edge of geometry is needed also to build the TDI observ- 
ables in such a way that laser phase noise and optical- 
bench noise are canceled effectively. In this section we go 
through all the elements of Fig. ^ an d discuss in detail 
how they are modeled in Synthetic LISA. In Sec. Ill Al wc 
describe the geometry of the LISA array, and the setup of 
the interferometric payload on each spacecraft; in Sees. 
Ill Bl and III CI we describe the response of the basic inter- 
ferometric observables to GWs and to the LISA funda- 
mental noises; last, in Sec. Ill Dl we give a rapid overview 
of TDI as used in LISA. 



A. LISA geometry and interferometry 

The motion of the LISA array is complex: at the qual- 
itative level, the three LISA spacecraft maintain a quasi- 
equilateral triangular configuration (where the arms stay 
equal to about 1%) trailing the Earth along its orbit in 
the plane of the ecliptic; at the same time, the constel- 
lation maintains an inclination of tt/2 — it /Q = n/3 with 
respect to the plane of the ecliptic (as measured from the 
normal of the instantaneous plane of the LISA constel- 
lation to the normal to the plane of the ecliptic), and 
it performs a cartwheeling motion, rotating around the 
normal to the instantaneous LISA plane with a rotation 
period of a year. This picture is realized in practice by 
placing the three spacecraft on eccentric, inclined solar 
orbits |2|. 

This pattern of motion improves the sensitivity of LISA 
to GW signals, making it more homogeneous over the 
sky (because the dependence of the antenna patterns to 
source position is averaged during the year), and improv- 
ing the estimation of source position and polarization 
(because the GW responses become modulated by the 
variation of the antenna patterns). This added sensitiv- 
ity comes at the price of complicating the GW response: 
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FIG. 2: Schematic diagram of the proof-mass and optical- 
bench assemblies within each LISA spacecraft (adapted from 
Ref. El). 

the modulations induced by the changing orientation of 
the LISA plane spread the power of originally monochro- 
matic GW signals, generating several sidebands at fre- 
quency multiples of 1/yr [lOj; furthermore, the relative 
motion of the detector with respect to the GW source 
introduces a time-dependent Doppler shift, which is the 
dominant effect for signals above 10 -3 Hz [the charac- 
teristic bandwidth of the Doppler shifting is ~ (QR/c)f, 
where / is the GW frequency and f2 = 27r/yr is the LISA 
orbital angular velocity]. 

When LISA is in operation, each spacecraft will ex- 
change laser beams with the other two, measuring the 
phase of the arriving laser beams with respect to the lo- 
cal lasers; the laser beams are bounced off freely-falling 
proof masses that are shielded by the spacecraft from 
most external disturbances, 1 so that they can serve as 
references for the measurement of GWs. To implement 
this measurement scheme, each spacecraft will carry 
two lasers, two proof masses, and two optical-readout 
schemes. FigureElpresents a schematic diagram (adapted 
from Ref. [ll|) of the proof-mass and optical-bench as- 
semblies within one of the LISA spacecraft, labeled "1" ; 
the other two spacecraft have identical setups. In short: 

1. the left-hand bench receives the laser beam from 
spacecraft 2, bounces it off its proof mass, and com- 
pares it with the local laser (without bouncing the 
latter) at the upper photodetector; 

2. via an optical fiber, the left-hand bench receives 
the right-hand-bench laser and compares it with 
the local laser (without bouncing the latter) at the 
lower photodetector; 



1 For instance, the spacecraft cannot completely shield the proof 
masses from cosmic rays. 
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FIG. 3: Schematic LISA configuration. The spacecraft are 
labeled 1, 2, and 3; each spacecraft contains two optical 
benches, denoted by 1, 1*, 2, 2*, 3*, as indicated. The 
unit vectors hi and light-path lengths Li connecting space- 
craft s and r are indexed by / = {1,2,3} for (s,r) = 
{(3, 2), (1,3), (2,1)}, and I = {-1,-2,-3} for (s,r) = 
{(2, 3), (3,1), (1,2)}. 

3. the left bench sends out the local laser (without 
bouncing it) to spacecraft 2, and (after bouncing 
it off the other side of its proof mass) to the right- 
hand bench. 

The operation of the right-hand bench (and indeed, of 
the benches on the other two spacecraft) is similar. [A re- 
cent candidate redesign of the optical benches ^3 would 
implement the comparison of the two lasers on the two 
benches of the same spacecraft by measuring their phases 
separately, doing away with the optical fiber, and then 
subtracting the measurements. For the purpose of ob- 
taining the laser-noise-free TDI signals (see Sec. IIIDf) 
this modification amounts only to a redefinition of the 
intra-spacecraft phase measurements 0, so m this pa- 
per, and indeed in Synthetic LISA, we refer to the older 
architecture.] 

In this setup, the physical observable of interest is the 
comparison of phase between the local laser and the in- 
coming laser, which carries information about the vari- 
ations induced by GWs in the inter-spacecraft optical 
path. The phase fluctuations of the lasers, however, are 
much larger than the GW-induced phase shifts, and must 
be subtracted before GWs can be resolved. In the last 
few years, a number of authors collaborated to develop a 
scheme ( Time-Delay Interferometry, or TDI) to subtract 
laser noise by carefully combining time-shifted series of 
the inter- and intra-spacecraft phase measurements; if 
the lasers are not phase-locked to a master (see the end 
of Sec. IIID|) . the intra-spacecraft phase measurements 
carry no information about GWs, but they do carry a 
combination of the phase noises from the lasers within 
each spacecraft. 

Because TDI has its origin in the techniques used to 
measure GWs by the Doppler tracking of distant space- 
craft 0, 0| , it prefers to describe the comparisons be- 
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tween laser beams in terms of fractional frequency dif- 
ferences rather than relative phase shifts (the two de- 
scriptions are exactly equivalent |lfij |. as they are re- 
lated by time integration). Thus, TDI represents the 
LISA readouts as basic Doppler observables: y s i r (t) is 
the fractional frequency difference at time t between the 
beam received at spacecraft r(eceiver) from spacecraft 
s(ender) and the local laser; and z s i r (t) is the analogous 
intra-spacecraft measurement on the same optical bench 
(thus, although it carries the index s it is in fact the 
fractional frequency difference between the two lasers of 
spacecraft r). In this paper, the index ^(ink) denotes the 
(oriented) LISA arm along which the laser was trans- 
mitted, according to the cyclical indexing I = {1,2,3} 
for (s,r) = {(3, 2), (1,3), (2,1)}, and I = {-1,-2,-3} 
for (s,r) = {(2, 3), (3,1), (1,2)} (thus, sgn(Z) = c.|,| r ). 
This spacecraft and link indexing is shown also in Fig. 01 
Note that our notation for the basic Doppler observables 
merges the two notations used in the scientific literature 
on first- generation TDI (yi r and z/ r ) and, more recently, 
on second-generation TDI (y sr and z sr ). TableQJshows a 
comparison (as it were, a Rosetta stone) of the notations 
used in various papers on TDI. In the next two sections 
we discuss the response of these basic Doppler observ- 
ables to GWs and to the noise sources present within 
each spacecraft. 



B. LISA response to gravitational waves 

In this section we give an expression for the GW re- 
sponse of the basic Doppler observables Vijif)- Work- 
ing in an inertial reference frame filled by a plane GW 
with propagation vector k and transverse-traceless grav- 
itational tensor h(x; t) = h(0, t — k ■ x) = h(i), we denote 
the positions of the three spacecraft by Pi(t). Following 
Estabrook and Wahlquist [T{|, we write the response of 
the inter-spacecraft LISA Doppler observable y s i r (t) to 
the plane wave as 

yfiTW = [! + &■ M*)] X (*i[*send ~ k-p s {t scnd )] 

-*l[t -*•&(*)]), (1) 

where i SC nd and p s {t sen d) are determined by the light- 
propagation equation t scnd = t - \p r (t) - p s (t scnd )\, 
where hi(t) is the oriented photon-propagation unit vec- 
tor hi(t) oc p r (t) — p s (t scn d), and where 



w W) = Mt)-h(ti)- ni (t) 

l[ 2(l-[fc-n i (<)] 2 )' 



(2) 



Equation gives the inter-spacecraft Doppler observ- 
able for laser-beam reception at time t on spacecraft r 
from spacecraft s, through link I. The two k ■ p products 
correspond to the retardation of the plane wavefront to 
the position of the two spacecraft, while the k-h products 
come m as geometrical projection factors M Equation 



JTJ is not singular for k = ±n;, because in that case the 
transverse-traceless tensor h is orthogonal to k and hi, 
and the *f?i go to zero. 

The light-propagation equation defines the effective 
armlength Li(t) experienced by light propagating from 
s to r, for reception at time t: 



Li(t) = \p r (t) - p s (t - Li(t))\. 



(3) 



Note that in general Li(t) ^ L_;(i). 

The response to GWs of the intra-spacecraft Doppler 
observable z s i r {i) is null, because the distance traveled 
by the intra-spacecraft beam is negligible for the GW 
amplitudes and wavelengths relevant to LISA. 2 



C. LISA response to fundamental noises 

In this section we give the response of the basic Doppler 
observables to the fundamental noise sources present 
within each spacecraft. Looking back to Fig. |2 we la- 
bel the left-hand and right-hand optical benches (and 
their lasers) as 1 and 1*, respectively (more generally, 
unstarred benches transmit into oriented arms with nega- 
tive indices) . Following Estabrook and colleagues i we 
denote the fractional frequency fluctuations of the laser 
on optical bench 1 as C%(t); these enter additively in 
the ?/23i measurement, together with the frequency noise 
from the laser on bench 2* of spacecraft 2, retarded to 
the time of emission: 



y%i sc (t) = CZ(t-L 3 (t))-C 1 (t) + 



(4) 



next, the 2/231 measurement is subject to noise due to fluc- 
tuations on the optical path of the beam incoming from 
spacecraft 2 (a combination of shot noise, pointing noise, 
and other optical-path noises), which we denote as 2/231! 
also, the velocity noise v% of the proof mass on optical 
bench 1 (i.e., its deviation from perfect free fall) induces 
an additional Doppler shift on the incoming beam (the 
local beam does not bring in any velocity noise, since it 
is not bounced on the local proof mass): 

2/23°™ = C^t-L,{t))-C 1 {t)+yll l (t)-2v 1 {t)-h z {t)+- ■ ■ 

(5) 

last, the random velocities V% and V\ of the emitting 
and receiving optical benches (which are several orders 
of magnitude greater than v\) induce additional Doppler 
shifts with the same temporal structure of laser frequency 



2/2 n 3 T(i) = C 2 *(t-L3W)-C 1 W+2/2 3 P iW-2^i(i)-«3W 
+ V 2 *(t-L 3 (t))-n 3 (t)-Vi(t)-h^(t). (6) 



2 If the lasers are not phase-locked to a master (see the end of Sec. 

ITTT51 . 
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reference 


geometry 


inter-spacecraft 
measurement 


intra-spacecraft 
measurement 

(on same bench of y 
with same indexes) 


link vectors 


armlengths 


this paper 


1 — > 2 3 cw 


y s l r (frequency) 


Z s [ r 


til (I > ccw, / < cw) 


L ; (/ > ccw, / < cw) 


AFT 1 000 r??l 


y ar (frequency) 




nm (along ccw) 


Lij (no distinction between 
ccw and cw) 


TAE 2000 [21], ETA 2000 [17], AET 

2001 [23], Hogan & Bender 2001 [21], 
Prince et al. 2002 [20], TDI whitepaper 

nnm nm act thai v°t A~\ 

2002 1 1 y J , Abl 2003 [24] 


z \l\r 


TEA 2002 [16], Tinto et al. 2003 [6] 


S|;i r (phase) 


T, flr (phase) 


Dhurandhar et al 2002 [25] 


U \ = >"3-21> ^2 =>l-23 
U 3 =^2-13' y l =^231 
^2 = "^312' y 3 = 


n/a 


ni« (along ccw) 


L|j (no distinction between 
ccw and cw) 


Cornish & Rubbo 2003 [8], Rubbo, 
Cornish, & Poujade 2004 [8] 




<J> Jr (phase) 




Z. (oriented) 


Hellings 2001 [18] 


1 — > 2 — ¥ 3 ccw 


y sr (phase) 


n/a 




y Jr (frequency) 


Krolak et al. 2004 [39] 


1 -> 2 -> 3 cw 


tim (along ccw) 


L (equal arms) 


Shaddock et al. 2003 a [4], 
TEA 2004" [5] 


i Jr (phase) 


t sr (phase) 


n/a 


L \l\ = L -b L \lV = L l 
(unprimed cw, primed ccw) 


Shaddock 2004 [26] 


L (oriented) 



TABLE I: A comparison of the phase-measurement and LISA geometry conventions used in the literature on TDI. In the cited 
references, A, E, and T refer to J. W. Armstrong, F. B. Estabrook, and M. Tinto. Notations are described with respect to the 
usage of this paper, with s = sending spacecraft, / = armlink, r = receiving spacecraft; "cw" and "ccw" refer to the progression 
of spacecraft or link indexes, as seen when looking at the LISA constellation from above (from ecliptic latitude 90° N); when 
indexes are shown in absolute values, only positive values are used. Tinto and Armstrong 1999 (not included in this table) has 
j/i = two-way ccw (a — > b — » a), ?/2 = two-way cw (a — > c ^ a). See www.vallis.org/tdi for updates to this table. 
a The semicolon ordered-delay notation was introduced in Shaddock et al. 2003 and TEA 2004. 
6 TEA 2004 uses to denote link vectors; it is ambiguous from the context whether these are ccw or cw. 



Along similar lines we derive the noise response of 
the intra-spacecraft measurement 23-21 on spacecraft 1, 
which contains the frequency noises from lasers 1 and 1* 
at time t, the random velocities of the optical bench 1 
and of its proof mass, and the frequency shift 771 upon 
transmission through the optical fiber (ultimately due to 
a component of the relative bench motions, V\ — V*): 

Z3-21 (*) = Ci(t) - Ct{t) + 2n 3 (t) ■ vi + 2h- 3 (t) ■V 1 +r n ; 

here we are ignoring time-delay effects along the fibers. 

Throughout the rest of this paper (and indeed, al- 
ways in Synthetic LISA) we take the optical-fiber noises 
and the optical-bench motions to be negligible. In fact, 
optical-fiber noise is removed in TDI by always using 
the z s ir observables in pairs such as (2:231 — ^3-2i)/2, 
(z 3 i 2 — zi_32)/2, and so on. One sees also that the optical- 
bench motions along the lines of sight (e.g., h-3 ■ Vi, 
h 2 ■ V*, and n 3 ■ V2) can be absorbed in the correspond- 
ing laser frequency noise variables (e.g., C\, CI, and C|), 
because they are appear in y^ sc and z" ; ° lsc with the 
same indices and the same evaluation times. Thus, if 
the TDI observables can successfully subtract laser fre- 
quency noise, they will also subtract the optical-bench 
motions, which are generally several orders of magnitude 



smaller. 

In writing Eqs. (Q, ©, and J7J), we have neglected 
also the offsets (up to several hundred MHz) between 
the center frequencies of the six LISA lasers, as well as 
the slow Doppler drifts resulting from the relative motion 
of the spacecraft (up to tens of MHz). In practice, the 
frequency offsets and Doppler drifts will be corrected by 
down-converting the photodetector output and tracking 
fringe rates using onboard ultrastable oscillators (USOs) 
D3 113 Although USOs introduce an important ad- 
ditional source of phase noise, their treatment is cumber- 
some, and we leave their modeling to a future version of 
Synthetic LISA. 

Under these assumptions, the simulation of the LISA 
noise response requires time series for 18 fundamental 
noise variables: the six proof-mass velocity noises along 
the line of sight (which we denote as pmi = n 3 ■ v±, 
pm 2 = hi ■ v 2 , pm 3 = hi ■ v 3 , and pm\ = n_ 2 • V\, 
ym\ = fi_3 • v^, pm 3 = h^i ■ v^), the six optical-path 
noises y°f r , and the six laser noises Ci and C* . (Note 
that our definition of the pmt differs by a factor — 1 from 
the definition used in Ref. [3.) 

The general expressions for jy";° lse and z" ; ° lse then be- 
come 
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and 



noises f C* s {t-L l (t))-C r (t)+yZ(t)-2pm r (t) if I > 0, 

It) — < (8) 

\a(<-ii(*))-c?(*)+Cr(*)- 2 i»n;(*) in <°> 



f C*(t)-C r (t) + 2pm*Jt) ifi>0, 



C r (t)-C;(t) + 2pm r (t) if Z < 0. 



We set standard levels for the 18 fundamental noises ac- 
cording to the noise budget discussed in the LISA pre- 
phase A report lj. Note however that Synthetic LISA 
allows all these prescriptions to be overridden. 

Laser Frequency Noise. We take each laser noise to 
be white, and to have a one-sided (square-root) spec- 
tral density of 30 Hz/VHz, which converts to a power 
spectrum of fractional frequency fluctuations by squar- 
ing and dividing by the square of the optical frequency 
~ c/(1064nm) = 2.82 x 10 14 Hz; thus, S% = 1.1 x 
10~ 26 Hz" 1 . We assume that the six laser noises are sta- 
tistically independent (the lasers need not be locked). 

Proof-Mass Noise. We take each proof mass to have 
white acceleration noise along the line of sight, with a 
one-sided (square-root) spectral density of 3 x 10 _ 15 m 
s -2 Hz -1 / 2 , which converts to a power spectrum of frac- 
tional frequency fluctuations pLTL Il9| by using the deriva- 
tive theorem for Fourier transforms, and dividing by c ; 
thus, S£ m = (3 x 10- 15 m s~ 2 H Z - 1 / 2 ) 2 /(4^ 2 / 2 c 2 ) = 
2.5 x 10 _48 [//Hz]- 2 Hz" 1 . We assume that the six proof- 
mass noises are statistically independent. 

Optical-Path Noise. We combine shot noise and 
beam-pointing noise on each optical bench into aggre- 
gate optical-path noises; we take these to be white dis- 
placement noises, with a one-sided (square-root) spec- 
tral density of 20 x 10~ 12 m Hz -1 / 2 , which converts to 
a power spectrum of fractional frequency fluctuations 
by using the derivative theorem for the Fourier trans- 
form, and dividing by c 2 ; thus, S° p = (20 x 10~ 12 m 
Hz- 1 '/ 2 ) 2 x (4tt 2 / 2 )/c 2 = 1.8 x 10~ 37 [f/Hz] 2 Hz" 1 . If 
the length of the LISA arms is different from the nominal 
value of 16.6782 s, we scale the optical-path rms noise by 
L;/(16.6782 s) to account for the 1/Lf power loss along 
the arms. 3 We assume that our aggregate optical-path 
noise enters the y s i r and z s i r observables in the same way 
as shot noise, and we further assume that the six optical- 
path noises are statistically independent. 



The variance of shot noise is inversely proportional to the number 
of photons received, which is proportional to the power received. 
Since power scales as l/L-f, rms shot noise must scale as L;. We 
assume that the remaining part of the aggregate optical-path 
noise scales in the same fashion. 




-yi-32(t)+ym(t+L2) y23i(t)+yi-n(t+L_ 3 ) 



FIG. 4: The four combinations of two basic Doppler observ- 
ables with emission or reception at spacecraft 1 and at time 
t. 



D. LISA TDI observables 
Time Delay Interferometry USEEHIIllEH 

mmiHElliililli^] isa technique t0 combine 

the basic Doppler variables y s i r and z s i r into compos- 
ite observables that are sensitive to GWs, but that are 
free of the otherwise overwhelming laser frequency noise 
(they are also free of optical-bench and fiber noise, as 
discussed above). To understand how TDI works, it 
is useful to tie the algebraic representation of the TDI 
observables to a visual picture of the path traveled by 
light between the LISA spacecraft. Looking at Fig. 0] 
let us consider the four combinations 2/231 (*) — 2/3—21 (i), 
J/231 (*)+ 2/123 (* + L 2 ), -?/i-32(£ + L-3) + 2/i23(£ + L2), and 
2/231 (t) + 2/1—32 (i + L_3); the two laser beams involved in 
each of these combinations are either arriving to space- 
craft 1 or leaving it at the time i; the retardations by the 
armlengths L2 and L_3 are needed because the Doppler 
observables are always labeled by the time of beam re- 
ception. Using Eq. JSJ, we see that for the first three 
combinations the contribution to laser frequency noise 
that is produced at time t is due to lasers 1 and 1*, and 
it is equal to —C\(t) + C*(t). In the fourth combina- 
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tion, 2/231 (t) + 2/1-32 {t + £-3), no laser noise is produced 
at time t, because the same laser is used as emitter and 
reference. For the first three combinations, the laser- 
noise contribution can be canceled by subtracting from 
the y s i r expressions given above the intra-spacecraft mea- 
surement (l/2)[z23i(t) — 23-21 (t)], whose laser-noise com- 
ponent is again —C\(t) + C*(t) (in fact, as noted above, 
each of 2231 (0 and — 23-21 (0 contains the combination 
— Ci(t) + C*(t), but the difference of the two 2 has the 
added advantage of canceling optical-fiber noise). 

Naturally, the laser noise that is produced at the times 
t— £3, t—L-2, and £+L_3 (in various combinations 



for the four y s i r expressions) is still not canceled. We 
see, however, that a combination of y s i r observables that 
corresponds graphically to a closed circuit would cancel 
laser noise completely; to build such a combination, wc 
need to delay the times of evaluation for the y s i r so that 
the tip or tail of each arrow meets another tip or tail 
(and only one!) at just the right time. A simple example, 
valid in the case when the Li are time-independent and 
all equal to L, traces a light path analogous to the path 
used in a Michelson interferometer (see the left panel of 

Pig- IB , 



[2/123 (* + £2) + y 3 -2i(t + L 2 + £_ 2 )] - [2/1-32^ + £-3) + 2/231 (t + £-3 + £3)] 

- 5(2231 (*) - Z3-21W) + |(^23l(* + £-3 + £3) - Z 3 _ 2 i(t + L 2 + £_ 2 )), (10) 

I 



where the two interfering light beams leave spacecraft 1 
at time t, and return at time t + 2L. The two double- z s i r 
subtraction terms are needed for the initial time of emis- 
sion of the two beams, and for the final time of arrival, 
while laser noise is self-canceling at the zero-angle cor- 
ners where the beams retrace their path, as mentioned 
above. Reordering Eq. pi)fl so that t is the final time of 
arrival of the beams at spacecraft 1, we get 

[2/3-21 (*) + 2/123, -2(f)] - [2/231 (*) + 1/1-32,3(0] 

- I 0*231,-33 (*) - ^3-21.2-2(0) + | 0*231 (0 - 23-2l(0): 

(11) 

where the comma notation y s i r ,d 1 d 2 ... denotes retardation 
by the armlengths L^, £d 2 j and so on. Laser-noise can- 
cellation works in this case because the length of the two 
paths 1 — > 3 — > 1 and 1 — > 2 — > 1 is the same (2L), so 
we can line up both the starting and the ending points 



of the two paths. If the arms (and hence the paths) were 
unequal, we would be left with residual laser noise origi- 
nating from the starting points of the two paths, as given 
by 

-(|Cl,-33(*) " 2-Cl,2-2(0) + (2^1*2-2(0 - 2-^1.-33(0) , 

The case of unequal (but constant) arms is tackled suc- 
cessfully by using new paths (1 — >3 — >1 — ^2 — >1 and 
1^2^1^3^1) each of which traces out both orig- 
inal paths (1 — > 2 — > 1 and 1 — > 3 — ► 1), but in opposite 
orders (see right panel of Fig. EJ. In this case, if we set 
the two paths to end at time t, the times of departure are 
both t — (L 2 + L-2) — (£-3 + £3), and the z s i r correction 
terms can cancel the noise emitted at that time, as well 
as time t. The corresponding TDI combination (known 
as unequal-arm-Michelson X, and first derived by Tinto 
and Armstrong is 



A = [2/1-32,32-2 + 2/231,2-2 + J/123,-2 + 2/3— 2l] _ [2/123,-2-33 + 2/3-21,-33 + 2/1-32,3 + 2/23l] 

— 2 (^3-21,2-2-33 — z 231, -332-2) — 2 ( z 3-21 — z 23l) (13) 
+ 5 (23-21,2-2 ~ z 231,2-2) + 3(23-21,-33 ~ 2 2 3i j _33) 

I 



where we omitted the dependence on / common to all the 
terms. 

Many TDI combinations are possible: all cancel laser 
noise, but each shows a different coupling to GWs and to 
the remaining system noises (known collectively as sec- 
ondary noises). As the understanding of TDI improved, 
the standard TDI observables evolved through various 



generations, capable of canceling laser noise for increas- 
ingly complex LISA geometries: 

First-generation TDI. Also known as TDI 1.0. The 

first-generation TDI observables [H SHI i3 cancel 
laser noise exactly in LISA configurations with unequal 
(but constant) arms, and Lk = £-&• Interferometric 
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Equal-arm Michelson Unequal-arm Michelson 




FIG. 5: Tracing the light paths in the Michelson and unequal- 
arm Michelson TDI combinations. 



combinations of various types are possible: 

The Sagnac-type observables (a, 0, 7) are sums of six 
basic Doppler observables, and they involve the difference 
between the Doppler shifts accumulated by light propa- 
gating around the LISA array in the two senses. Thus, 
the Sagnac-type observables use all the LISA laser links 
in both directions. A fully symmetric Sagnac observ- 
able (£) is considerably less sensitive than most others 
to GWs with frequencies at the lower end of the LISA 
band; it was suggested [2lJ that the comparison between 
the power observed in £ and in the other TDI variables 
could be used to discern a stochastic GW background 
from instrumental noise. The observables built from six 
Doppler variables are also known as six-pulse combina- 
tions, because their response to an impulsive plane GW 
consists of six separate pulses. 

Eight-pulse combinations involve sums and differences 
of the Doppler shifts measured along four of the six LISA 
laser links. The unequal-arm Michelson observables {X, 
Y, Z) use both links of two arms; as discussed above, 
they can be interpreted as measuring the phase differ- 
ence accumulated by light traveling (twice, in opposite 
orders) along the two arms of a Michelson interferometer 
centered in one of the spacecraft. Perhaps for this reason, 
and in analogy with ground-based GW interferometers, 
a single unequal-arm Michelson observable (generally X) 
is often used in LISA data analysis to compute expected 
detection rates and parameter-estimation accuracies. 

More eight-pulse combinations can be formed: the bea- 
con observables (P, Q, R) use only the two links depart- 
ing from one of the spacecraft, and both links along the 
opposite arm; the monitor observables (E, F, G) use only 
the two links arriving at one of the spacecraft, and both 
links along the opposite arm; last, the relay observables 
(U, V, W) use one departing link and the adjacent arriv- 
ing link at one of the spacecraft, together with both links 
along the opposite arm. The eight-pulse combinations 
can be considered as LISA contingency modes, because 
they are available even if one or two of the laser links fail. 
Note however that all six lasers must still be available to 
build the intra-spacecraft observables z s i r required for 
the eight-pulse combinations, except in the case of the 



unequal-arm Michelson observables: one of these can al- 
ways be built even if one or both lasers directed along 
one of the arms happen to fail. 

Dhurandhar and colleagues j2f| proved that the space 
of all the first-generation TDI observables can be con- 
structed by combining four generators, which they iden- 
tify in a, P, 7, and (. Prince and colleagues [2(j showed 
how to diagonalize the cross noise spectrum of the gen- 
erators to obtain three observables (^4, E, and T) with 
uncorrelated noises. The three optimal observables A, E, 
and T are written as sums and differences of a, (3, and 7, 
and when used in combination they achieve the optimal 
S/N for GW sources at any frequency in the LISA band. 

Modified TDI. Also known as TDI 1.5. Shaddock [H 
recently pointed out that the rotation of the LISA array 
introduces a difference in the armlengths experienced by 
beams traveling in the corotating and counterrotating 
directions (i.e., Lk 7^ £-&)• Furthermore, this difference 
becomes much larger if we take into account also the 
orbital motion of the array around the Sun [5j. Some 
of the first-generation observables (the A"-type, P-type, 
i?-type, and [/-type combinations), cancel laser noise 
also for Lk 7^ if time delays for the appropriate 

oriented arms are used [as we have already arranged, 
for instance, in Eq. i|14fl ]; these observables can be in- 
terpreted as tracing light paths that enclose vanishing 
areas. Conversely, the first-generation observables that 
trace light paths that enclose a finite area (such as a, 
0jl, and C) are equivalent to Sagnac interferometers 
|28|. and must necessarily be sensitive to the rotation 
of the array, which shows up as a spurious phase dif- 
ference between the lasers, originating from the starting 
points of the light paths. The Sagnac observables can be 
modified by means of a finite-difference procedure anal- 
ogous to the change undergone between the equal-arm 
and unequal-arm Michelson combinations (see Fig. [SJ) , so 
that the modified Sagnac observables have null enclosed 
area, and cancel laser noise 0, The resulting com- 
binations [ai, and 0:3, which generalize a, /3, and 7; 
and Ci, C2, and £3 p|, which nonuniquely generalize (] in- 
clude twice as many y s i r variables as the first-generation 
combinations (i.e., they are 12-pulse observables). 

Second-generation TDI. Also known as TDI 2.0. The 

motion of the LISA array introduces not only a direc- 
tional dependence of the armlengths, but also a time de- 
pendence, as first recognized by Cornish and Hellings 0. 
In this case, the order of the TDI retardations becomes 
important: for instance, if the armlengths are constant, 
then 

t t 2— 2 = t — L-2 ~ L2 = t — L2 — L-2 = t,-22 (14) 

but if they are not (as signaled by a semicolon index 
notation), then 

ta-2 = (t- L_ a (t) - L 2 (t - L- 2 )) 

^(t-L 2 (t)-L_ 2 (t-L 2 ))=t ; -22. (15) 
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More generally, the semicolon notation represents the re- 
tardation chain rule 



t-, dl ...d n = t - L dn (t) - L An _ x (t - L dn {t)) 
- L dn _ 2 (t - L dn (t) - L dn _ t (t - L dn (*))) - 



(16) 



where the rightmost retardation index is applied first, 
using the armlength L dn (t) ; the next-to-rightmost retar- 
dation index is applied second, using the partially re- 
tarded armlength L dn l (t — L dri (t)), and so on. Taylor- 
expanding the armlengths, and retaining only the zeroth- 
order and first-order terms, we get 



t- dl ... dn — t — L dn — L dn _ 1 — L dn _ 1 L d 



L d „_„ (L„ 



L 



(17) 



where for ease of notation we have dropped the (t) de- 
pendence common to all the armlengths. As discussed 
in Refs. 0, 0, H|, the eight-pulse TDI observables can 
be generalized, once again by a procedure akin to finite 
differentiation, to 16-pulse observables that cancel laser 
noise up to first order in the Taylor-expanded armlengths; 
for the LISA orbital parameters, this is enough to cancel 
laser noise to a level below the secondary noises. Accord- 
ing to the notation of Ref. , Xi , X2 , and A3 generalize 
X, Y, and Z; Pi, P 2 , and P3 generalize P, Q, and R; 
Ei, E 2 , and E 3 generalize E, F, and G; and U\, U 2 , and 
U3 generalize U, V, and W. The Xf. observables can be 
interpreted as expressing the difference in laser phase be- 
tween beams propagating along two paths whose Taylor- 
expanded total lengths differ only by terms proportional 
to Lk or to higher derivatives; 4 the residual laser noise is 
then a sum of expressions similar to 

C k . A (t) - C k , B (t) C k {t) x [t. A - t. B ] 

~ Cfe(t) x 0[L and higher derivatives] . (18) 

As for the Sagnac-type observables, the 12-pulse modified 
observables ai, a 2 , ocs, (1, C2, and (3 can already cancel 
laser noise to a level below the LISA secondary noises: 
the residual laser noise is of order L and higher, but the 
specific combination of L k involved turns out to be small 
for the LISA orbit. 

Although historically the TDI observables were de- 
rived by combining time-shifted combinations of the ba- 
sic (one-way) Doppler measurements y s ir{t), they can 
also be written as combinations of one-way and two-way 
Doppler measurements, generated by locking five of the 
six lasers to the remaining one, as described by Tinto 



and colleagues |6j ; the resulting expressions contain fewer 
terms, are still noise-canceling, and have the same re- 
sponse to GWs. 



III. IMPLEMENTATION AND USAGE OF 
SYNTHETIC LISA 

Synthetic LISA is an object-oriented CH — h library built 
to mirror the idealized structure of Fig. ^ each block 
in the figure corresponds to one or more CH — I- classes 
|29j |. which implement its functionality. The Synthetic 
LISA workflow follows this object-oriented structure, fa- 
cilitating targeted investigations that compare multiple 
configurations of one object (for instance, one of the fun- 
damental noises, or the GW source), while all others are 
kept fixed. Here is an example of a typical Synthetic 
LISA session. 

1. Create an instance of a LISA geometry (LISA) class 
with the desired orbital parameters. 

The LISA classes provide the geometrical quanti- 
ties Pi(t), hi(t), and Li(t) needed to assemble the 
LISA GW and noise responses described in Sees. 
Ill Bl and III CI They account for the aberration ef- 
fects caused by the finite speed of light and by the 
spacecraft motion intervening between the events 
of pulse emission and reception [Eq. 

There are different levels of complexity at which the 
motion of the LISA array, discussed in Sec. Ill AI 
can be modeled in a simulation of the LISA sci- 
ence process; correspondingly, increasingly sophis- 
ticated TDI observables are needed to cancel laser 
noise once the added complexity is taken into ac- 
count. In Synthetic LISA, these levels correspond 
to different derived classes 5 of the base class LISA. 
The simplest such class, OriginalLISA, models a 
stationary, nonorbiting constellation, used implic- 
itly in the development of first-generation TDI. The 
most realistic, Eccentriclnclined, models the ec- 
centric orbits of the spacecraft up to second order 
in the eccentricity (see App.[SJ; the resulting time 
dependence of the armlengths creates the ne- 
cessity of second-generation TDI for effective laser- 
noise suppression. 

2. Create instances of a LISA noise class (Noise) for 
the 18 fundamental-noise time series defined in Sec. 
1/761 tuning noise parameters if so desired. 

SyntheticLISA can generate pseudorandom noise 
sequences that approximate closely the standard 
laser, proof-mass, and optical-path noises specified 
in Sec. Ill CI alternatively, the package can import 



4 The finite differencing procedure adopts the compound paths A 
= I+II and B = II+I, where the paths I and II must contain the 
same links, in different orders; then tin ~ till — I X II — 11x1 = 

(Ei LOiEj L n .) - e, imXEj Li t ) = 0.' 



J In CH — h, a derived class inherits the data content and behavior 
of its base class, and can add enhancements or customizations. 
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the noises as sampled time series, which might have 
been generated with other tools, or even measured 
experimentally. The treatment of the LISA noise 
processes is crucial to the interpretation of Syn- 
thetic LISA simulations, and is discussed in de- 
tail in App. [5] ln short, the representation of 
noise is adequate if the noise-generation Nyquist 
frequency /& is set comfortably higher than the 
highest frequency at which one wishes to analyze 
the TDI noise responses, but of course lower than 
the Nyquist frequency used to sample the TDI ob- 
servables, to avoid aliasing. 

3. Create an instance of a GW source class (Wave) of 
the desired type and parameters. 

The Wave classes provide the GW polarization com- 
ponents h+(t) and h x (t), which are assembled into 
the transverse-traceless metric perturbation h(t) 
according to the polarization convention described 
in App. ^ Synthetic LISA contains simple Wave 
classes (such as SimpleBinary for monochromatic 
binaries), which can be modified easily to yield 
more complicated signals; the package can also im- 
port h x (t) and h x (t) as sampled time series. 

4. Create an instance of a LISA TDI class (TDI), feed- 
ing it the LISA geometry, LISA noises, and GW 
source objects previously created. 

The base class TDI defines a complete set of first- 
generation, modified, and second-generation TDI 
observables, according to the expressions of Refs. 
[Hill for first-generation TDI, and of Refs. HQ 
for modified and second-generation TDI. 6 The de- 
rived classes TDIsignal and TDInoise implement, 
respectively, the LISA response to GWs [y^ r from 
Eq. QJ], an d to the fundamental noises [y";° lse and 
z s°r S ° lrom E( l s - © an d ©]■ Users can easily de- 
fine additional TDI observables, using Table [Q to 
rewrite the expressions in the literature in terms of 
the y s i r and z s i r Synthetic LISA observables. 

5. Last, use the TDI objects to generate a time se- 
ries of the TDI observables and write it to disk or 
memory. 

No C++ programming and compilation is needed to use 
Synthetic LISA, since the functionality of the package 
can be accessed very easily from the scripting language 
Python |30| . either interactively, or with short scripts. In 
fact, the Synthetic LISA session described above would 
translate to a handful of lines in Python. Refer to 
the Synthetic LISA manual |3l| for detailed information 
about the usage and implementation of the package. The 



The primed link indices of Refs. U Q correspond to positive 
indices in this paper. 



manual documents also the successful validation of Syn- 
thetic LISA's output against analytical expressions of the 
TDI observables for both noise and signals. 

IV. NUMERICAL EXPERIMENTS WITH 
SYNTHETIC LISA 

We now present the main scientific results of this pa- 
per: an investigation of laser phase noise suppression 
for flexing LISA array orbits with first- and second- 
generation TDI [Sec. lIV A*| : an analysis of the armlength- 
determination accuracies required for effective laser-noise 
suppression [Sec. UVTj] : and an estimation of quantiza- 
tion and telemetry bitdepth needed for the phase mea- 
surements y s i r and z s i r . While significant as they stand, 
these studies are meant also to exemplify the kind of 
system-characterization inquiries that becomes possible 
with advanced LISA simulators. 

Except where otherwise specified, all the power spectra 
displayed in this section were computed as periodograms, 
reducing spectral leakage and fluctuations by dividing 
one-year-long time series into partially overlapping seg- 
ments (in number of either 1024 or 2048, depending on 
the specific test), triangle- windowing each segment, and 
averaging the resulting power spectra (see, e.g., Ref. |32fl). 
Thus, all the spectra of this section represent average ef- 
fects: slightly different requirements on laser-noise power, 
armlcngth determination, and phase-measurement quan- 
tization might be needed to achieve the same suppression 
performance homogeneously across the year. 

A. On the necessity of second-generation TDI 

As recognized by Cornish and Hellings |3j, the eccen- 
tric and inclined orbital motion of the LISA spacecraft 
introduces a time variation in the armlengths of order 
10~ 8 s/s [sec Eq. I|A3|) of Apd.IA"]; as a consequence, the 
first-generation and modified TDI observables fail to can- 
cel laser frequency noise completely. Using the graphical 
interpretation of TDI given in Sec. Ill Dl we would see 
that the interferometric circuits synthesized by the ob- 
servables fail to close exactly. The laser-noise residuals 
arise from the starting points of the paths, and they are 
of the form 

Sd = |[^ J (0-C* 7 (t)]-|[C i;J (t)-C' i;J (t)] (19) 
~ \[C*{t) + Ci(t)](t.j - tj) = i[C*(t) + CifflSt, 

where / and J denote time-ordered path retardation 
chains. Using the Fourier derivative theorem and assum- 
ing white, uncorrelated laser noises, we get 

|<5C< (/)| 2 /l^| 2 = 2vr 2 / 2 ^ 2 , (20) 

for frequencies up to the laser- noise bandlimit. For the 
modified TDI X observable, 5t ~ 10~ 6 s, so laser noise 
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FIG. 6: Imperfect cancellation of noise with modified TDI 
X for flexing LISA array (Eccentriclnclined, £o — rjo = 0). 
The top curve plots the perfect-cancellation noise target, ob- 
tained by setting the laser noise to zero; the intermediate and 
bottom curves show the noise spectra resulting from imper- 
fect laser-noise cancellation for nominal and 0.3 x rms laser 
noise. A curve with O.lx rms laser noise would be essentially 
indistinguishable from the perfect-cancellation target, and is 
not plotted here. The spectra are computed from one year's 
worth of X data sampled at 1 Hz (with 1-s noise-generation 
timestep), averaging over 2048 data subsegments. 
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FIG. 7: Reduction in amplitude S/N due to imperfect laser- 
noise cancellation for modified TDI X, with realistic LISA 
spacecraft orbits (Eccentriclnclined, £o = ??o = 0) and with 
nominal laser noise, 0.3x rms laser noise, and O.lx rms laser 
noise. Because of the linear frequency axis used in this graph, 
the bump between 1 mHz and 10 mHz (which tops at about 
1.3) appears unimpressive, but it spans a scientifically impor- 
tant frequency range. (From the same data as Fig. [(J) 



is canceled by less than 160 c!B for / > 2 mHz; for the 
second-generation TDI X\ observable, St ~ 10~ 10 s, so 
laser noise is canceled comfortably by more than 160 dB 
throughout the LISA band of good sensitivity. In this 
section we discuss the results of Synthetic LISA simu- 
lations carried out to investigate and substantiate these 
analytic arguments. 

Figure shows the spectrum of secondary noise plus 
residual laser noise (top curve) versus the spectrum of 
secondary noise alone for the modified TDI X observ- 
able (bottom curve), computed using realistic eccentric 
and inclined LISA spacecraft orbits; the excess noise 
is evident between 1 mHz and 10 mHz, and within 
the noise nulls at the frequency multiples of 1/(2L). 
The intermediate curve shows the effect of reducing the 
laser noise to 0.3 times its nominal rms spectral density 
1.1 x 10~ 26 Hz . A separate simulation was performed 
by reducing laser noise to 0.1 times its nominal value: the 
resulting spectrum is essentially indistinguishable from 
the secondary-noise-only curve. 

The reader might be puzzled by the flatness of the 
first-generation TDI noise curves at low frequencies, as 
compared to the f~ 2 dependence of proof-mass noise and 
of the often-seen LISA sensitivity curves. The flatness 
is caused by the time-delay structure of first-generation 
TDI observables, which contain, as it were, a finite- 
difference time derivative, with a low-frequency power 
transfer function proportional to f 2 . On the other hand, 
the sensitivity curves plot a ratio of (rms) noise response 



to GW-signal response, with the latter decreasing as 
/ 2 at low f req uency for first-generation TDI observables 
such as X 17]. 

Figure shows the reduction caused by residual laser 
noise in the (amplitude) S/N for monochromatic sources, 
computed as the square-root ratio of the imperfect- 
cancellation and secondary-noise-only spectra. The loss 
of sensitivity appears significant (up to ~ 30%) between 
1 mHz and 10 mHz, and even more so around the 1/(2L) 
harmonics. However, an improvement in laser noise sta- 
bility by a factor of about three would be sufficient to 
erase the S/N-reduction bump at lower frequencies, and 
to shrink considerably the S /N-reduction peaks at higher 
frequencies. An improvement by a factor of ten would es- 
sentially eliminate the need for second-generation TDI, 
as estimated analytically in Ref. Q . 

By contrast, Fig.[S]shows that essentially perfect laser- 
noise cancellation is achieved with the second-generation 
TDI observable X\, with residual laser noise (bottom 
curve) several orders of magnitude below the secondary 
noises. For the Sagnac observable Q\ (which, strictly 
speaking, belongs to the set of modified TDI observ- 
ables) . laser noise is still canceled by more than one order 
of magnitude, except at the first d null near 6 x 10~ 2 
Hz. 
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FIG. 8: Cancellation of laser noise with second-generation 
TDI Xi and modified TDI fi for flexing LISA array 
(Eccentriclnclined, £o = T)o — 0). For Xi, the noise spec- 
trum obtained by including nominal laser noise matches ex- 
actly the perfect-cancellation noise target obtained by setting 
the laser noise to zero (top curve); the bottom curve, several 
orders of magnitude below, shows the spectrum of residual 
laser noise. For fi, residual laser noise sits only one order of 
magnitude below the perfect-cancellation noise target. The 
spectra are computed from one year's worth of Xi and £i 
data sampled at 1 Hz (with 1-s noise-generation timestep), 
averaging over 2048 data subsegments. 



B. On the required accuracy of armlength 
determination 

We now turn to simulating the laser-noise residuals 
resulting from the implementation of TDI using an im- 
perfect knowledge of the armlengths, and therefore of the 
TDI delays. In this case, residuals are created at all the 
delayed times that appear in the TDI observables, and 
not just at the starting point of the interferometric cir- 
cuits. Graphically, the reason is that the tail of each y s i r 
arrow (determined by a physical light-travel delay) does 
not precisely match the head of the preceding arrow (de- 
termined by a nominal TDI delay affected by armlength- 
measurement error). At each such point, the residual has 
the form 



6C s ~C s (t)6Li(t). 



(21) 



The total residual noise is a somewhat complicated func- 
tion of the TDI observable under consideration. 

Tinto and colleagues |(| find that an armlength ac- 
curacy of ~ 30 m (~ 100 ns) would be needed for ef- 
fective laser-noise subtraction with first-generation TDI. 
They also estimate how often the armlength measure- 
ments would have to be updated, by computing the 
timescale for the time-dependent armlengths to change 
by an amount equal to the required accuracy; for realistic 
LISA orbits, this timescale varies substantially through 



the year, but it can be as low as 10 s. 

In the course of the LISA mission, armlengths might 
be determined by means of orbital-dynamics models that 
are periodically updated by ranging measurements, ei- 
ther between the spacecraft, or to Earth. It was recently 
suggested j^H that the TDI observables do not need to 
be assembled in real time aboard the spacecraft, but that 
they can be obtained in postprocessing from time series of 
the y s i r and z s i r measurements sampled at limited rates 
(~ 1 Hz) and telemetered to Earth. If that is the case, 
the accuracy of ranging is probably a secondary issue, 
since even poor measurements can be fitted a posteriori 
to very accurate models of the LISA orbits. In fact, it 
was recently proposed that the ranging information 
can be obtained directly from the y s i r and z s i r measure- 
ments, by minimizing the integrated noise power in the 
TDI observables as a function of the orbital parameters 
of the LISA spacecraft. Because of these reasons, the 
problem of determining the accuracy required for ranging 
measurements is not well defined in the context of post- 
processed TDI. In this section we concentrate instead on 
the accuracy required for the real-time onboard compu- 
tation of the TDI observables. 

The simplest real-time treatment of the TDI delays 
consists simply of keeping the armlengths fixed to their 
last measured values, which are updated at time intervals 
T. The resulting requirements on the ranging measure- 
ments are rather constraining: for modified TDI X, our 
simulations show that marginally acceptable laser-noise 
cancellation is obtained with measurements repeated ev- 
ery 8 s with 2-m (rms) accuracies (assuming independent 
errors). Indeed, the piecewise-constant armlength model 
does a very poor job of following the dominant linear 
time dependence of the armlengths. 

A better treatment, which requires very little sophisti- 
cation in the onboard logic, consists of extrapolating lin- 
early from the latest two armlength measurements, which 
are again repeated at intervals T. The left panel of Fig. 
[5]shows that, for modified TDI X, 100-m (rms) accurate 
measurements, repeated only every 4096 s, yield residual 
laser-noise suppression to better than a factor of six be- 
low the case of perfect armlength knowledge (where some 
residual laser noise is present because of the LISA array 
flexing; see Fig.|SJ). Every successive n-fold improvement 
in the accuracies yields an n 2 -fold improvement in laser- 
noise suppression. 

Remarkably, taking ranging measurements more of- 
ten has the effect of worsening laser-noise suppression at 
low frequencies. To understand why, consider Eq. I|21|l . 
which implies that in the Fourier domain the laser-noise 
residual is given by the convolution of the laser-noise 
derivative with the armlength error. The rapid repeti- 
tion of measurements introduces high-frequency power 
in the armlength-error time series with a typical band- 
width of 1/(2T), which then causes the leakage of power 
from high frequencies (where C s is much larger) to the 
low-frequency end of the LISA spectrum. This behavior 
can be observed in the left panel of Fig. by comparing 




FIG. 9: Imperfect cancellation of laser noise with modified TDI X (left panel) and second-generation TDI Xi (right panel) 
due to imperfect knowledge of the armlengths in a flexing LISA array (Eccentriclnclined, £o = T]o = 0). The topmost 
curves show the result of using perfect armlengths: thus, the X\ curve shows secondary noise only, while the X curve shows 
secondary noise plus the residual laser noise due to using modified TDI observables with a flexing LISA array. All other curves 
show the level of residual laser noise for linearly-extrapolated armlengths (see main text) with different single-measurement 
errors AL rms and intervals T. Error-laden armlength measurements are simulated by adding a Gaussian-distributed, zero-mean 
independent deviates to the correct values of the six Li. The low- frequency flattening of the laser-noise residuals is caused 
by power leakage from high frequencies when the bandwidth of the armlength-error time series is comparable with the LISA 
measurement bandwidth. The spectra are computed from one year's worth of X and Xi data sampled at 2 Hz (with 0.5-s 
noise-generation timestep), averaging over 1024 data subsegments. 
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FIG. 10: Imperfect cancellation of laser noise with second- 
generation TDI X\ due to quantization of the phase measure- 
ments using riquant bits (see main text). The strongly sloping 
curve plots only the secondary noises, while all other curves 
show the level of residual laser noise with different quantiza- 
tion depths (the numbers shown do not include the additional 
1 + 3 bits needed for the sign and to avoid saturation). The 
spectra are computed from one year's worth of X and X\ 
data sampled at 2 Hz (with 0.5-s noise-generation timestep), 
averaging over 1024 data subsegments. 



the laser-noise residual curves corresponding to measure- 
ments repeated every 4096 s and every 64 s. By contrast, 
the maximum acceptable spacing of the measurements 
is set by the timescale for relevant quadratic changes 
in the armlengths: for a typical armlength acceleration 
a ~ (27r/yr) x 10~ 8 s/s = 2 x 10~ 15 s _1 , the time required 
to accrete an error ~ 100 m is ~ y/2 x lOOm/a = 18,000 
s. 

The right panel of Fig. EH shows that the armlength 
accuracy requirement for the second-generation TDI ob- 
servable X\ is not substantially different, with 100-s 
(rms) accuracy achieving laser-noise suppression by a 
(power) factor of about five, and successive n-fold accu- 
racy improvements yielding n 2 -fold suppression improve- 
ments. However, considerably better accuracy is needed 
if laser noise is to be canceled also within the X\ nulls 
at 1/(2L) and multiple frequencies. The leakage effect 
discussed above is more important in the case of second- 
generation TDI, where at low frequencies secondary noise 
declines as a positive power of /, and can intersect the 
leakage plateau if measurements are not taken sparsely 
enough. 

C. On the quantization of phase measurements 

Our last numerical experiment in this paper is con- 
cerned with estimating the number of effective bits that 
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must be obtained and recorded for the phase measure- 
ments y s i r and z s \ r , and then either transmitted between 
the spacecraft to perform TDI in real time, or transmit- 
ted to Earth to perform it in post processing. 

Similar, less extensive experiments have been per- 
formed by J. W. Armstrong The underlying physical 
problem is that laser noise must be represented faithfully 
enough to allow its cancellation by several orders of mag- 
nitude. Thus, we expect the spectral characteristics of 
laser noise, such as its bandwidth at the output of the 
phasemeter, and its magnitude relative to the secondary 
noises, to play into the answer to our question. Presum- 
ably, considerable telemetry bandwidth can be saved by 
whitening phase noise prior to transmission, in such a 
way that the quantity of (Fourier-space) information rel- 
ative to secondary noise is approximately constant at all 
frequencies. For the purpose of our estimates, we adopt 
the crude whitening scheme implicit in dealing directly 
with fractional-frequency fluctuations; in this paper we 
assume laser noise to be white for these. 

We quantize phase measurements by dividing each 
y s i r and z s \ r (before assembling the TDI observables) by 
a fiducial fractional-frequency-fluctuation level given by 
the nominal rms value of laser noise (i.e., 1.05 x 10~ 13 , 
assuming noise bandlimited at 1 Hz), truncating the re- 
sulting values to n quant bits to the right of the binary 
point, and then multiplying again by the fiducial level. 
The actual counting of bits must include one additional 
sign bit, and a few bits to the left of the binary point (we 
take three, which is adequate to make the truncation of 
Gaussian-distributed noise statistically insignicant). Fig- 
ure 1101 shows the results of our simulations for second- 
generation TDI X\. an n quant between 32 and 34 (and 
hence a total number of bits between 36 and 38) is needed 
to lower the level of residual laser noise resulting from 
quantization to a level comfortably below the secondary 
noises in the LISA measurement band. The requirement 
is less strict (ri quant between 28 and 30) for modified TDI 
X. 

More definitive simulations of the effects of measure- 
ment quantization should include less idealized models of 
phase noise at the output of the phasemeter. Note also 
that the simulations presented here do not address the 
interplay between quantization and the implementation 
of fractional- filtering interpolation, used in postprocessed 
TDI [33] to approximate the values of y s i r and z s i r at the 
TDI delays between recorded samples. 

V. CONCLUSIONS 

We have described three numerical experiments on the 
implementation of TDI in LISA, which were performed 
with Synthetic LISA, a simulation of the LISA science 
process that can generate synthetic time series of fun- 
damental noises and GW signals, as they appear in the 
laser-noise-canceling TDI observables. Our conclusions 
were presented in brief in SecQJ and described in detail 
in Sec. II VI We have also discussed the theoretical model 



that underlies Synthetic LISA and provided details of its 
implementation, as needed to understand the results of 
our numerical experiments. 

The structure and programming style used for Syn- 
thetic LISA allows for vast extensions and improvements. 
Among others, we plan to include explicitly the addi- 
tional time series required for calibration of the onboard 
ultrastable oscillators ^(|, and to model explicitly the 
measurement errors at the photodetectors. We are in 
the process of making Synthetic LISA available |36j as 
a public-domain software package, to foster the involve- 
ment of the wider GW community in research on the 
interface between scientific goals and technical require- 
ments for LISA, on the tradeoffs and improvements that 
can be made in the implementation and operation of the 
mission, and on the development of novel analysis tech- 
niques for the LISA data. In the spirit of open-source 
design, we expect the LISA and GW communities to pro- 
vide their own useful additions to Synthetic LISA, such as 
more realistic models of the noises and of the spacecraft 
subsystems, and additional GW source modules. For this 
purpose, we have designed Synthetic LISA as a modular 
and easily extensible C++ package, with a user-friendly 
Python frontend for easy scripting and prototyping. 

The investigations that can be carried out with state- 
of-the-art simulators such as Synthetic LISA include: 

Performance characterization and architecture 
trade-off studies. Synthetic time series supplement 
analytical results in the allocation of subsystem noise 
budgets and in the determination of the final sensitivity 
for specific GW sources, providing a high-level analysis 
tool for system engineering, and helping the formulation 
of technical requirements from the desired LISA science 
goals. 

An example was the recent study [32l of detection 
prospects for the GW signals from compact stellar ob- 
jects inspiraling into the supermassive black holes at 
the centers of galaxies, with the purpose of determining 
whether the LISA noise floor would need to be lowered 
to guarantee a minimum number of such detections. For 
this study, time series for h + and h x were produced us- 
ing the Glampedakis-Hughes-Kennefick quasiadiabatic 
orbit integrator (3^, and then fed to Synthetic LISA, 
which computed the corresponding time series of TDI 
observables; these were used to derive the expected S/Ns 
for the capture sources. 

Noise analysis and vetos. Synthetic time series can 
be used to study real-LISA features of the instrumen- 
tal noises, such as nonstationarity, noise increments due 
to faulty subsystems, or (perhaps most important) the 
level of cancellation of laser phase noise by TDI under 
different LISA geometries, armlength-measurement tol- 
erances, and other TDI characteristics. The numerical 
experiments presented in this paper represent a first step 
in the numerical validation of TDI as implemented for 
LISA; more detailed studies will undoubtedly become 
necessary as additional details about the actual imple- 
mentation of TDI become available. 
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FIG. 11: Conventional definition of the GW polarizations + 
(dashed) and x (solid) for various ecliptic latitudes j3 and 
longitudes A. Figure excerpted from Ref. |3?j |. 
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APPENDIX A: 



GEOMETRIC CONVENTIONS IN 
SYNTHETIC LISA 



A recent example was the use of Synthetic LISA |2J| 
to validate a new approach to the determination of the 
LISA armlengths, whereby the noise power in the TDI 
observables is minimized as a function of the armlengths. 

Development of data-analysis algorithms. The syn- 
thetic time series produced by this simulation have con- 
sistent signal structure and noise correlations across all 
the TDI combinations. Thus they can be used to test 
algorithms for use on the real LISA data, such as the 
separation of stochastic GW backgrounds from LISA in- 
strumental noises |2lJ , the matched-filtering detection of 
quasiperiodic signals [39| . and so on. Synthetic LISA 
provides a streamlined module to filter GWs through the 
LISA TDI response, allowing easy interfacing to existing 
GW data or GW-modeling applications. GW data an- 
alysts using Synthetic LISA to generate simulated LISA 
data will also be able to exploit the library of GW signals 
being assembled at the Mock LISA data archive [40||. 



In this Appendix we give a brief overview of the Syn- 
thetic LISA conventions for the LISA orbits and for the 
geometry of GW sources, with the purpose of facilitating 
the comparison between the synthetic signals generated 
by Synthetic LISA, and those created with other methods 
(and in particular with the LISA Simulator |8(). As dis- 
cussed in Sec. Ill Al the motion of the LISA array baricen- 
ter is approximately contained in the plane of the ecliptic; 
for this and other reasons |5j, Synthetic LISA employs a 
Solar-system-baricentric ecliptic coordinate system (SSB 
frame), setting the x axis toward the vernal point. 

The most accurate description of the LISA mo- 
tion available in Synthetic LISA (the LISA class 
Eccentriclnclined) models the eccentric orbits of the 
spacecraft up to second order in the eccentricity e. For 
these orbits, the dominant (and evolving) differences be- 
tween the armlengths are caused by the flexing motion of 
the array 0] due to orbital eccentricities. Following Ref. 
the SSB coordinates pi of spacecraft i are given by 
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where a, = 3tt/2 - 2{i - 1)tt/3 and e = 0.00964838, 
yielding an effective L ~ 16.6782 s. These spacecraft 
orbits are mapped to those used in the LISA Simulator 
Q by setting 770 = ft, Co — 37r/2 — k+A, where n and A are 
the parameters defined below Eqs. (56) and (57) of Ref. 
0, and by choosing sw < in the Eccentriclnclined 
constructor, which has the effect of exchanging spacecraft 
2 and 3. 

The armlengths experienced by light propagating 



along the arms can be found by solving Eq. For 
efficiency, Synthetic LISA employs the lowest-order ap- 
proximation 
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where Si = {£o,£o + 47r/3,^o + 27r/3}. The amplitude 
of the flexing correction is about 7.5 x 10~ 2 s, or 0.5% 
of the nominal LISA armlength; the rate of change of 
the armlengths is about 1.5 x 10~ 8 s/s, which requires 
second-generation TDI to yield sufficient cancellation of 

laser phase noise. 

All the Synthetic LISA GW source objects (Wave) 

share the same geom etrical setup, which follows the con- 
ventions of Ref. .'!!) . At the position x in the SSB frame, 
the spatial part of the transverse-traceless metric per- 
turbation associated with a plane GW can be written as 

h(t) = h+(t - k ■ x)e + + h x (t - k ■ x)e x ; (A4) 
here the functions h + (t) and h x (t) express the two po- 



larization components of the wave at time t, measured at 
the origin of the SSB frame. For a GW source at ecliptic 
latitude j3 and ecliptic longitude A, the unit propagation 
vector k is 

k= — (cos (3 cos A, cos (3 sin A, sin 0). (A5) 

The two polarization tensors e+ and e x that appear in 
Eq. l|A4jl are defined without loss of generality as 

/ 1 \ /0 1 0\ 

e+ = E • -1 • E T , e x = E ■ 1 • E T , 

\ o o o / \oooy 

(A6) 

where the orthogonal matrix E, 



sin A cos ip — cos A sin j3 sin ip — sin A sin ip — cos A sin (3 cos ip — cos A cos f3 
— cos A cos ip — sin A sin f3 sin ip cos A sin ip — sin A sin (3 cos ip — sin A cos [3 



(A7) 



cos 13 sin ip 



cos j3 cos tp 



sin j3 



expresses an Euler rotation sequence, whereby the (3 and 
A terms can be understood as enforcing the transversal- 
ity of the GW, while the polarization angle tp encodes a 
rotation around the direction of wave propagation, — k, 
setting the convention used to define the two polariza- 
tions. The polarizations corresponding to ip = are 
shown in Fig. 1111 for various source positions in the sky. 
The positional parameters (3, A, and ip are mapped to 
the parameters 9, <p, and ip used in the LISA Simulator 

[1] by setting f3 = ir/2 - 9, A = </>, and ip = -ip. 

The standard monochromatic-binary Wave object, 
SimpleBinary implements the GW signal 



-(*) 
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.4 



(1 + COS 2 l) X C0S(27T ft + <po) 

(2 cos i) x sin(27r/£ + </>o) 



(A8) 



where A is the common amplitude, t is the inclination 
angle, / is the GW frequency observed in the SSB frame, 
and (po is the phase at t = 0. The standard value of 
A is (2mim2/ ! dR) with mi, W2 the two masses, d the 
luminosity distance, and R the orbital separation (the 
common amplitude ho used in Ref. |39| differs by a fac- 
tor of two, ho — 2 A, absorbed in h^ and h+). We 
have found excellent agreement (see, e.g., Fig. IT2)l be- 
tween the time series of TDI observables derived from 
SimpleBinary and the output of the LISA Simulator, v. 
2.0 [8j using Newtonian, c. 



deviates, 7 which are then interpreted as the sampled val- 
ues at times t n = nAt (for n = 0, 1, . . .) of a continuous 
random process. The process is assumed to be bandlim- 
ited below .ft, = 1/(2 At): by the sampling theorem (see, 
e.g., Ref. [33I ), the value of the noise can then be recon- 
structed exactly at any intervening time t by convolving 
the sampled sequence with the interpolating kernel 



sinc[7r(i — t n )/At] 



sin [tt (i - t n )/At] 

[TT(t-t n )/At] 



(Bl) 



Since the sine kernel has infinite time extent, it must 
be replaced in practice by an approximated interpolation 
scheme that involves a finite number of samples. A vast 
class of such schemes, including the linear and polynomial 
interpolators implemented in Synthetic LISA, can be for- 
mulated as the convolution of the sampled sequence with 
an interpolating kernel that is (in some sense) an approx- 
imation to the sine. 

The tradeoff in the approximation is between the num- 
ber of samples used to interpolate and the sharpness of 
the spectral response. The correct sampling of a ban- 
dlimited process preserves all the spectral information 
below the Nyquist frequency, but it populates Fourier 
space with infinitely many replicas of the original spec- 
trum, centered at frequencies k/ At (for k — ±1,±2, . . .). 
The effect of sine interpolation is to multiply this compos- 
ite spectrum by the Fourier transform of the sine, which 



APPENDIX B: SYNTHETIC NOISE IN 
SYNTHETIC LISA 

In Synthetic LISA, pseudorandom white noise is cre- 
ated by generating a sequence of uncorrelated Gaussian 



7 Independent, uniformly distributed dev iate s are obtained from 
Luescher's lagged Fibonacci ge nerator |41| . as implemented in 
the GNU Scientific Library 14 j : the Box-Muller transform 14:1 
is then used to convert the uniform deviates to Gaussian deviates. 
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S(f ) Monochromatic binary at/= 1.944 mHz, lst-gen. X 
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FIG. 12: Comparison between the outputs of the LISA Sim- 
ulator, v. 2.0 [3 using Newtonian, c and of Synthetic LISA for 
the response of first-generation TDI X to a monochromatic 
binary with / = 1.944 mHz, j3 = \ = ip = Q, t~ 7r/2, and 
amplitude appropriate for a (0.5 + O.O33)M0 binary at 0.1 
kpc. We show the nonaveraged spectra of the unwindowed 
signals in a neighborhood of the central frequency. To com- 
pare the fractional-frequency-fluctuation output of Synthetic 
LISA to the nominal-strain output of the LISA Simulator, we 
have multiplied the strain spectrum by the square of the nom- 
inal armlength (10 10 m) to get a displacement spectrum; con- 
verted to a velocity spectrum using the derivative theorem for 
Fourier transforms; and converted to a fractional-frequency- 
fluctuation spectrum by dividing by c 2 . 



S(f) [a.u ] Spectral noise distortion due to interpolation 
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FIG. 13: Noise distortion for 0.5 Hz-bandlimited white noise, 
with different interpolation schemes. Higher-order interpola- 
tion corresponds to a sharper transition at the bandlimit fre- 
quency and to lower ripples and deeper valleys between them 
(in this graph, the valleys are cut off by spectral leakage from 
the main platform at power = 1). 



is a perfect square window of height 1 and width 1 / At, 
centered at / = 0. Thus, sine interpolation achieves per- 
fect signal reconstruction by selecting only the original 
spectrum and deleting all unwanted replicas. Practical 
schemes with kernels of finite extent cannot have such a 
sharp frequency response, so they distort (i.e., amplify 
or suppress, depending on frequency) the original spec- 
tral content in the passband below /;,, and they allow 
some of the power of the unwanted spectral replicas to 
creep back into the interpolated process (either directly, 
if the process is sampled with a sufficiently high Nyquist 
frequency, or by aliasing to frequencies i n th e passband) . 

These effects can be observed in Fig. which shows 
a spectrum of pseudorandom white noise, generated with 
a timestep of 1 s, and resampled to a timestep of 0.1 s, 
using no interpolation (i.e., defaulting to the nearest 1-s 
sample), using linear interpolation, and using Lagrange- 
polynomial interpolation of order 4, 8, and 32. In all 
cases, power begins to drop before the nominal bandlimit 
frequency of 0.5 Hz, but the drop is sharper and closer to 
0.5 Hz for higher-order interpolation methods. Spurious 
power above the bandlimit frequency appears as ripples 
between the /& harmonics: the height of the ripples de- 
creases with the interpolation order, while the valleys 
among the ripples become wider. In Fig. 1131 the valleys 
appear to be cut off by a common downgrading enve- 
lope; this is an artifact of spectral estimation, due to the 
residual leakage from the platform below the passband; 
spectral leakage also smears out to a finite height the 
nulls at the fb harmonics. 

In Synthetic LISA, interpolated pseudorandom white 
noise is used to stand in for the standard laser phase noise 
of Sec. Ill CI The standard proof-mass and optical-path 
noises, which have colored spectra, are approximated by 
applying simple digital time-domain filters to the uncor- 
related deviates, before interpolation. Namely, the finite- 
difference filter y[nAt] — x[nAt] — x[(n — l)Ai] (with x 
the original noise sequence) has power transfer function 
|1— exp(27ri/A<)| 2 = 4sin 2 (7r/Ai), and is used to approx- 
imate the standard f~ 2 proof-mass noise. The damped- 
integrator filter y[nAt) — ay[(n — l)Ai] + a;[nAt] (with 
a = 0.9999, to control the DC component of y) has power- 
transfer function ~ (1/4) sin _2 (7r/Ai), and is used to ap- 
proximate the standard f 2 optical-path noise. 

The resulting pseudorandom noises have power spec- 
tra that adhere very faithfully to the nominal curves, 
except at frequencies comparable to ft, where the effect 
of interpolation is that noise power is not cut off sharply, 
but rather drops off smoothly (if rapidly) , with nulls at 
the fb harmonics. For the optical-path and proof-mass 
noises, the effect of interpolation is compounded by the 
effect of the finite-difference and finite-integration time- 
domain filters, whose transfer function near fb is pro- 
portional to sin [nf /(2fb)] rather than / ±2 . We con- 
clude that the pseudorandom noises can be accurate rep- 
resentations of the standard LISA noises of Sec. Ill Cl and 
therefore can be used to study the noise response of the 
TDI observables, as long as we take into account the ef- 
fects of interpolation and filtering at frequencies compa- 
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rable to /;,. Because TDI is essentially a linear operation, 
the results at lower frequencies will not be affected. Us- 
ing linear interpolation (the Synthetic LISA default), it 
is probably safe to draw conclusions from the TDI results 
at frequencies < /&/5; using higher-order interpolation, 
it becomes possible to push inferences to higher frequen- 
cies. 

This discussion of filtering and interpolation applies 



also to noise objects provided by the user as sampled 
time series, as long as the sampled noise can be consid- 
ered bandlimited below its nominal Nyquist frequency. 
See Ref. j3^| for a related discussion of the use of inter- 
polation in reconstructing the TDI observables on Earth 
from the y s i r and z s i r data, sampled onboard at a limited 
rate that can be transmitted affordably to Earth. 
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